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GLOBAL EXISTENCE AND ASYMPTOTIC STABILITY OF 
SMOOTH SOLUTIONS TO A FLUID DYNAMICS MODEL OF 
BIOFILMS IN ONE SPACE DIMENSION 

R. BIANCHINI 0 & R. NATALINI* * 


Abstract. In this paper, we present an analytical study, in the one space 
dimensional case, of the fluid dynamics system proposed in [4] to model the 
formation of biofilms. After showing the hyperbolicity of the system, we show 
that, in a open neighborhood of the physical parameters, the system is totally 
dissipative near its unique non vanishing equilibrium point. Using this property, 
we are able to prove existence and uniqueness of global smooth solutions to the 
Cauchy problem on the whole line for small perturbations of this equilibrium 
point and the solutions are shown to converge exponentially in time at the 
equilibrium state. 


Keywords. Fluid dynamics models, dissipative hyperbolic equations, biofilms, global 
existence, asymptotic stability. 

1. Introduction 

A biofilm is a complex gel-like aggregation of microorganisms like bacteria, algae, 
protozoa and fungi. They stick together, attach to a surface and embed themselves 
in a self-produced extracellular matrix of polymeric substances, called EPS. 

In this paper, we study a fluid dynamics model, introduced in [3], to describe the 
space-time growth of biofilms. This model was built in the framework of mixture 
theory, see [T2j or [T], and conserves the finite speed of propagation of the fronts. For 
simplicity reasons, the model describes a biofilm in which there is just one species of 
microorganisms, or better, all species are lumped together, but it can be extended 
to other situations. It has been derived starting from the equations for mass and 
momentum conservation, and some physical constraints and assumptions about the 
behavior of the biological aggregates and their interaction with the surrounding liquid. 
Here we assume that the complex structure of biofilms is described by four different 
phases: bacteria B(x,t), extracellular matrix EPS E(x,t), dead cells D(x,t) and 
a liquid phase L(x,t). The quantities B,E,D,L are the volume fraction of each 
component, then B,E,D,L £ [0,1]. Since we are dealing with a one dimensional 
model, then we have x £ M, t > 0. For simplicity reasons, we assume that B , E, D 
have the same transport velocity vs- The reaction terms are indicated by r<j>, with 
$ = B,E,D,L. Imposing the total balance of mass and momentum for each phase 
<h, we can write the model, see [4] for all details. This model was originally proposed 


° Dipartimento di Matematica, Universita degli Studi di Roma ’’Tor Vergata”, via della Ricerca 
Scientilica 1, 1-00133 Rome, Italy. 

* Istituto per le Applicazioni del Calcolo ”M. Picone”, Consiglio Nazionale delle Ricerche, via dei 
Taurini 19, 1-00181 Rome, Italy. 


1 



2 


R. BIANCHINI 0 & R. NATALINI* 


in all space dimensions, and in the present case of one space dimension, is given by 
a system of six partial differential equations, which read: 

[ d t B + d x (Bv s ) = T B , 

3 t E + d x (Evs) = Te, 
d t D + d x (Dv s ) = T D , 

(1.1) d t L + d x (Lv L ) = T l , 

<9 t ((l - L)v s ) + <9 X ((1 - L)vg ) = -(1 - L)d x P - 7^(1 - L) 

+ (M - T l )v l - Mvs ; 

_ d t (Lv L ) + d x {Lvl) = -L8 X P - (M - T l )v l - Mv s . 

To reformulate our model in a more suitable form, we assume the following volume 
constraint: 


( 1 . 2 ) L = l-(B + E + D), 

that is the assumption that the mixture is saturated and no empty space is left. In 
addition to the balance mass of each component, we also have the total conservation 
of the mass of the mixture by the following assumption: 

(1.3) T# + Tg + T^ + IY = 0. 

The mass constraint in m states that the mixture is closed, i.e. there is no net 
production of mass for the mixture. According to [5], the reaction terms are given 


by: 


(1.4) 

Vf> = K b BL - Kfj Ik 

(1.5) 

r E) = OlKe>B -- KjyD 

( 1 . 6 ) 

T e = E ePE — cE. 


The birth of new cells at a point depends on the quantity of liquid available in 
the neighborhood of the point, that is why the birth term in Tg is a product 
between the volume fraction B of active cells and the volume fraction L of liquid. In 
this way, the mass production term V B is the difference between a birth term and a 
death term, where the second is proportional to the fraction B of bacteria, with 
rate k B . The death term in the expression of T B gives rise to a creation term of the 
mass exchange rate for dead cells T B , with a proportional coefficient cc, since a part 
of the active cells goes into liquid when the cell dies. In T^, we also find a natural 
decay of dead cells with a constant decay rate kjy. The EPS is produced by active 
cells in presence of liquid and therefore the production term will be proportional to 
BL, where kE is the growth rate of EPS. There is also a natural decay of EPS with 
rate e. To conclude this explanation about the mass exchange terms, we choose T l 
in order to enforce condition (El. See again [5] for more details. 


Let us now further simplify the equations of system ED- First, by the volume 
constraint, the equation for the fraction L is no more necessary. Summing the 
equations for B, E, D , L in ED and using the volume constraint El, we obtain 

(1.7) <9 X ((1 - L)v s + Lv l ) = 0 , 
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which establishes that the derivative of the average of the hydrodynamic velocity 
vanishes. Since we are in one space dimension and since all phases B , E, D , L and 
velocities vs,vl vanish at infinity, by using equation (11.71) we can express the unknown 
liquid phase velocity in function of the solid phase velocity, and then we can eliminate 
the equation for the liquid phase from system (EH). Precisely, we have: 

, n , L - 1 

(1-8) v L = vs- 

Summing up the fifth and the sixth equation of system EH and using equality (11.81) , 
we have an explicit expression of the spatial derivative of the pressure term P, which 
solves 


(1.9) d x P = —7^(1 — L) — d x ((l - L)vg + Lvl). 

Then, substituting equation (11.81) in ESD, we have 


( 1 . 10 ) 


d x P = -7<9x(1 -L)- d x 


1 -L 


L 


-v s 


After some calculations, we obtain an equation for the solid phase velocity vs, and, 
finally, system EH becomes a system of only four equations, which are 

d t B + d x (Bv) = T b , 
d t E + d x (Ev) = r E , 

(1.11) { d t D + d x (Dv) =T d , 


d t v + d x 


(3 L - 2)v 2 
2 L 


+ 7 (L + ln( 1 - L )) 


r l-m ^ 

L{l-L ) V ~ Vv ’ 


where, from now on, v := vs- Obviously, we still have the biological constraints EH) 
and E3D, which close system (11.111) . 

In this paper, we aim to give a first analytical result to the Cauchy problem for 
this system on the real line. The paper is organized as follows. In Section 2, we 
determine the region of hyperbolic symmetrizability of system (11.111) . that ensures the 
local existence of smooth solutions (see 0, EJ). Then, we establish some conditions 
in order to have a totally dissipative source. Finally, in the last section, we use this 
dissipative property to prove that, if we have initial data in a small neighborhood - in 
// 2 -norm - of the unique equilibrium point u = (B,E,D,v), these smooth solutions 
are global in time. In the last part of the paper, we analyze the asymptotic behavior of 
our solutions by using energy estimates. More precisely, we show that these solutions 
decay exponentially, in the P 2 -norm, to the equilibrium point u. 

The proof of global existence of smooth solutions for system (11.111) takes inspiration 
from [5J, but here we do not make use of any strictly convex entropy. As we will 
discuss later, system (11.111) is ’’totally dissipative”, which means that the source term 
has a particular dissipative property. Our proof of global existence for the solutions 
of (11.111) is based on this totally dissipative property, which means that dissipation 
prevents, for small initial data, the formation of shocks and ensures global existence 
of the solutions. 


2. The quasilinear hyperbolic dissipative system 


Let us rewrite now system (11.111) in a compact form, in order to study its hyperbolic 
symmetrizability, which provides the local existence of smooth solutions (see El, 0)- 
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Let us set 


u = 


/ B \ 
E 
D 

\ v 


The flux function for system (11.111) is 

( 


( 2 . 1 ) 


and 


F{ u) = 


Bv 

Ev 

Dv 


\ 


(3 L - 2)v- 
2L 


+ 7 (L + ln(l - L)) j 


( 2 . 2 ) 


G{ u) = 


( Tb\ 

r £ 
r D 

V r «/ 

is the vector of the reaction terms. With this notation, we can write system (11.111) 
in the compact conservative form 

d x u + d x F( u) = G( u). 

To study the properties of (11.111) . we need to compute the Jacobian matrix 

\ 

n 11 n fi 

(2.3) A(n) = 


f V 

0 

0 

B 

0 

V 

0 

E 

0 

0 

V 

D 

V v 

V 

V 

(3L-2) 

L 


of the flux function F( u). We can now write (11.111) as a quasilinear system of partial 
differential equation, in the following form 


(2.4) 


<9tu + J 4(u)9 x u = G(u). 


2.1. Hyperbolicity and symmetrizability. The main hypothesis of the theorem 
on local existence of smooth solutions for a general hyperbolic quasilinear system 
in form m is symmetrizability (see i, 0), namely the existence of a symmetric 
positive definite matrix Alo; such that AoA is also symmetric. 

Assume L / 0 and L / 1, and set 

. L7 i ' 2 

(2.5) r, : = 


v 

w 


(1 -L) 

It easy to see that, if 77 > 0, then a positive symmetrizer for the model in (11.111) is 
given by the matrix 


A 0 (u) = 


( V 

0 

0 

0 ^ 

0 

Bv 

E 

0 

0 

0 

0 

By 

D 

0 

V 0 

0 

0 

B ) 


( 2 . 6 ) 
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Remark 2.1. By the existence of a symmetrizer, system {2.4V , 12.3\) . \2.2\) is 
symmetrizable hyperbolic in the following domain 


(2.7) 


W = < u = ( B , E, D, v) e[0,lfxl:i)> 0, 0 < L < 1 


Let us find now the domain of simple hyperbolicity of our model & m, m- 
We calculate the eigenvalues of A(u) in \2.3\) . which are the following: 

(1) Ai, 2 = v; 

(2) A 2 ,3 = (2L ~ 1)t; ± \/(l — L)A; 
where 


( 2 . 8 ) 


. L'y v 2 

A ^-i >0 ’ 


in order to have real eigenvalues. Moreover, observing \2.5\ ) we note that 


A > 7] > 0. 


This inequality is completely in accordance with the general theory about hyperbolic 
symmetrizable system, because (strong) hyperbolicity follows directly from hyperbolic 
symmetrizability. 

Therefore, always if rj > 0, we can write the explicit expression of domain of 
symmetrizability (and then hyperbolicity) for system {2.4V , {2.3\) . {2.2)) . which is 
(2.9) 


W = { u = (B, E, D, v) e [0, l] 3 x R : - 


L 3 / 2 7 V 2 
(1 -L)V2 


< v < 


L 3 / 2 7 V 2 
(1 -L)V2’ 


0 < L < 1 


Finally, let £1 a convex open subset of W. This set Ll will be our domain of work 
for system {2-4V , 112.3)) . {2.21) . 


Let us set now Al(u) = -Ao-A(u). We can write system (12.41) . (12.31) . (12.21) in the 
following symmetric form 


(2.10) A)(u)<9 t u + A^vifdxU = A)G(u), 

with Aq,A, G given respectively in (12.61) . (12.31) . (12.21) . 

From symmetrizability, we have local in time smooth solution in H 2 (W) for system 
m , dm m with initial data in the same space, thanks to the theorem of 
local well-posedness of the Cauchy problem for quasilinear hyperbolic symmetrizable 
systems (0). 


2.2. Total dissipativity of the source. Here we are going to show that, under 
some assumptions, the source term is totally dissipative and, then, under some 
smallness hypothesis on the initial data uo, these solutions are also global in time. 
This dissipative property, in fact, allows to estimate the Ed 2 - norm of solution u and 
to prove that this norm is bounded for all t > 0, without the use of strictly convex 
and strictly dissipative entropy, as done in [5]). 

Let us give now our definition of Hyperbolic Symmetrizable Totally Dissipative 
System. In order to clearly write this definition, we consider a general one-dimensional 
nxn hyperbolic symmetrizable system in the compact form (12.41) . with u£SlC M ra , 
where H is a convex open set and Al(u), G(u) are smooth enough in fL 
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Definition 2.1 ((D)-Condition). Assume that system 112.101) is hyperbolic 
symmetrizable in a set W. Let Ll be a convex subset ofW, which contains a unique 
equilibrium point u for 12.1(A) , i.e. G{u ) = 0. We say that system \2.10\ ) is totally 
dissipative in Q, if there exists a matrix D = D(u,u) £ M nxn such that, for every 
u £ Ll, we have 

• G(u) = D(u,u)(u — u); 

• AqD(u, u) is negative definite. 


Let us check the totally dissipative property of system (11.111) . First of all, we 
determine the coordinates of the point in which the reaction terms vanish. Actually, 
we have only one isolated equilibrium point u, which is inside of the domain of 
hyperbolicity. 

To see that, set 

b = ('i_M/('i + £*£ + *£*£ Y 

y Kb J / \ kn ek B J 

Therefore 


( 2 . 11 ) 


( 1 \ 

k B kD 


( 1 \ 

k B k B 


u = B 


ek B 

ak B 


\ 


kN 

0 




a 


i , ak D 
k N 


kpks 

ek B 


V 


ek B 

ak B 

kN 

0 




with a := kB k ^ D • It has been assumed that B,E,D,L £ (0,1). Then, from (12.111) . 
in order to have positive volume fractions B , E, D, we have to assume 


(2.12) k B > k D . 

Let us observe now that inequality (12.91) . which describes the region of hyperbolic 
symmetrizability of (11.111) . is verified near the equilibrium point u. Then, if we take 
u in a small neighborhood of this point, we are within the domain of hyperbolic 
symmetrizability. For this reason, to show our results of global existence of smooth 
solutions, the idea is to consider an initial data uo in a convex and compact 
neighborhood of this equilibrium point u, then we set 


(2.13) 


f l = B r (u), 


for a fixed r > 0. Let us write the matrix D, as in (D)-Condition, for the source 
term of system (11.111) . Then 





G{ u) = 

D(u,u)(u- 

u )> 

with 


( 

—Bk B 

—Bk B 

—Bk B 




eE(L-B) 

-e(BL+BE) 

eE 

(2.14) 

5 

C 

Cl 

II 


k R Jf) 

BL 

0 

kpjB 




B 

B 



V 

0 

0 

0 


° \ 
0 

0 

r l-m 

L(l—L) / 


In order to show that system in (12.101) . (12.61) . (12.31) . (12.21) verifies the (D)-Condition 
in a small neighborhood of u, we have only to prove that matrix AqD(u, u) is negative 
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definite. Using the expression of u, we rewrite the expressions of the matrices 


A 


o(u) = 


kp'y 

(k B —k D ) 

0 

0 

0 


0 

k E (k B —kp) 

0 

0 


0 

0 

fcjy 7 

a(k B -k D ) 

0 


0 \ 
0 
0 
B ) 


and 


D{ u,u) = 


( ~Bk B 
k E (L-B) 
ak B 


V 


0 


—Bkp 
-e — Bk E 
0 

0 


—Bk E 

—Bk E 

—kjy 

0 


0 

0 

0 

fc|M 


kp{k B —kp) ) 


Using the Routh-Hurwitz conditions ([lT)j) on - A ° D ') + ( A ° D ) ; we a i m to establish 

the sign of AqD(u, u). We have 
(2.15) 




—Bksko 

!f[e(L-B)-k D B} 
!f[k N -Bk B \ 


¥ [e(L -B)- k D B] ^ i k N - Bk B \ 


ek B {e+k E B ) 

k E _ 

ekftB 

2 

o 


eksB 

2 

k 2 

— Hn. 


0 
0 
0 

k 2 B MB 


. -P _ 

kD(k B —k D ) / 


Obviously, we are interested to the internal 3x3 matrix which contains the main 
diagonal of the originale 4x4 matrix. Its characteristic polynomial is P( A), and its 
general form is 


(2.16) 


P( A) — A 3 + aiA' + a 2 A + 03 . 


Here, the coefficients (i = 1, 2, 3) are all real. We require the conditions on the cq, 
such that the zeros of P(A) have Re A < 0. The necessary and sufficient conditions 
for this to hold are the Routh-Hurwitz conditions [IOJ. For the cubic equation, A 3 + 
aiA 2 + a 2 A + <23 = 0, the Routh-Hurwitz conditions for Re A < 0 are given by the 
following [RH] conditions: 


(2.17) 


[RH] < 


a\ > 0 ; 
a 3 > 0; 

0102 — a 3 > 0 . 
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In our case, the coefficients of (12.161) are the following: 
(2.18) 

e 2 k B ^ k% 


a\ = Bk B k B + 


. + eB + 

a 


N . 


e 2 Bk B k B eBk B k 2 N ek B k B k 2 N Bk B k B (eBk B + ek B + k B kN + e 2 
a2 = -;-1-1-h 


kE Oi OL 

B 2 k 2 B (k 2 D + e 2 ) k 2 D (k 2 N + e 2 )' 


e 2 Bk B k B kisr | e 2 Bk 2 B k B k‘^ l | eBk B k B k 2 N {Bk B + e + k B ) 

° 3 “ 2 k^ + ofa + 2 a 

eBk B k B kN(e + Bk B ) eBk B k B {tBk B kN + eBk B k B + k B k 2 N ) 

H 4 4 

e 2 k B k 2 D {k 2 N + B 2 k B ) k 2 N {e 2 k 2 D + e 2 B 2 k B + B 2 k^k B ) 

. 4/^e 4a 

Therefore, we can state our total dissipation result. 

Proposition 2.1. Assuming that f2.12\) holds. Then, if [RH] condition is valid for 
oiifl 2 j «3 (E 18\) . system \2.KA) . \2.b\) . 112. A) . \2.2\) satisfies (D)-Condition, and 
so it is totally dissipative in a neighborhood of its equilibrium point. 

Now, the first condition of (12.171) for a± in (12.181) is always verified. Regarding the 
second two conditions, it can be verified that they hold for some e, a, k B ,k B , kE, k n, 
in particular for the values of the parameters in the following table, which was used 
in [ 1 ]: 


Table 1. A list of (dimensional) parameters |T]. 


Par am. 

Value 

Unit of meas. 

Indications 

k B 

r> 

1 

O 

t—H 

00 

1 /sec 

Bact. growth rate 

kE 

12 • 10 ~ b 

1 /sec 

EPS growth rate 

k B 

2•10 Y 

1 /sec 

Bact. death rate 

kN 

1 ■ 10 -d 

1 /sec 

Dead cells consumption 

e 

1.25 ■ 10 -Y 

1 /sec 

EPS death rate 

a 

0.25 

dimensionless 

coeff. liquid dead-cells 


More generally, if we consider a costant a E [0.5,1.5], we can restrict our attention 
at the class of coefficients such that 

(2.19) e = 1.25 • 10 - 7 ~ 10~ 7 , kj\r = 10ae, ,kE = lOOae, k B = 2ae, k B = 70ae. 

By this reduction, condition for 03 (the third condition of (12.171) ) becomes a second 
degree polynomial inequality 

-3.366175 • 10 10 a 2 + 4.1829869 • 10 10 a - 8.303370950 • 10 9 > 0, 

that can be easily solved, and the inequality holds for a E [0.5,1]. Moreover, by 
a simple numerical verification, it can be seen that the second condition in (12.171) 
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is also true in this interval. Then, if we take the coefficients as in (|2.19l) . then 
(D)-Condition holds ture for a £ [0.5,1]. 

3. Global existence of smooth solutions for small initial data and 

ASYMPTOTIC BEHAVIOR 

Let us present now the proof of the global existence of smooth solutions to system 
(11.111) with initial data in a small neighborhood, in H 2 - norm, of the equilibrium 
point. We want to follow the proof of global existence of smooth solutions for weakly 
dissipative hyperbolic systems with a convex entropy in [5] and |2j. Using (Di- 
Condition, we can write the L 2 -norm estimates of the local solution u and of its 
first and second derivatives. Then, global existence of smooth solutions follows from 
these estimates. We underline that, in order to get these estimates, we use only 
the dissipative property of the source term - (D)-Condition - and we do not use 
any strictly convex and strictly dissipative entropy for system (II. Ill) , which is an 
assumption of the theorem of global existence in [5j. 

Now, if we take the initial data in a small neighborhood of the equilibrium point u, 
and we assume that our parameters are such that the (D)-Condition is verified, 
we are able to prove global existence of smooth solutions to model (11.111) . Moreover, 
using the totally dissipative property of (D)-Condition, we can prove that this 
global solution will decay to the unique equilibrium point. In the following, let us 
state our global existence and asymptotic behavior result. 

Theorem 3.1. Consider system and its unique equilibrium point u in \2.11\ ) 

and assume that \2.12\ ) holds. If this system satisfies (D)-Condition, then there 
exists a 0 < 5 < r such that, if \ \uq — u\ [2 < 6, then there is a unique global solution 
u with initial data uq, which verifies 

u-u£C{[0,+oo),H 2 {R))nC 1 ([0,+oo),H\R)) 

and 

(3.1) sup 0 <t<+oo ||(u- u)(t)\\l + [ \\(u-u)(t)\\l dr < C(5)\\u 0 - u\\ 2 , 

Jo 

where C(5) is a positive constant. Moreover, the global solution u decays exponentially 
in H 2 -norm to the equilibrium point u, i.e. 

(3-2) \\u- m||h 2 (m) < C\e *||uo - m||h 2 (R) , 

where C\,/3 are positive constants. 

In accordance with the hypothesis, we work in a neighborhood U of the equilibrium 
point u. Let us introduce the new unknown 

(3.3) w := u — u. 

Then, the equation in (I2.1UI) becomes 

(3.4) A)( w + u)dtw + j4i(w + u)<9 x w = j4oG(w + u). 

To prove global existence for (11.111) . we follow the approach proposed in [11], see 
also [6j and [5], and we introduce the functional 
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Proposition 3.1. Let T > 0, and assume that there exists a smooth solution w of 
12.1(A) . 12.61) . 112. 3)) . 112.2)) in [0,T], Then, there exists e > 0 and C > 0 such that, if 
N 2 (T)<e, 


(3.6) iVf(T) < C(Nl( 0) + 1V 2 3 (T)). 

In a classical way, the first part of Theorem \3.1\ about global existence and 
uniqueness of smooth solutions follows directly from Proposition \3.1\ (see m, mi). 
The following Lemma is the first step of the proof of global existence. 

Lemma 3.1. If N 2 (T) < e, then 

(3.7) Nq(T) < Ci (IVf(0) + 1V 2 3 (T)). 

Usually, to state an estimate as (3ZD, a function of convex entropy is used, but 
here, in our proof of Lemma 13.11 we do not use anything but the dissipative property 
of system (11.1111 . 

Proof. Using (D)-Condition and (GOD, let us consider system (11.1111 in the following 
symmetric form: 

(3.8) A)( w + u)<9t\v + Ai(w + u)c? x w = Ao(w + u)D(u, u)w. 

In the previous equation, the new reaction term is 

(3.9) L>i(u, u, w) := A 0 (w + u)L>(u, u). 

Therefore, we have 

(3.10) A)( w + u)5 f w + Al(w + u)5 x w = L>i(u, u, w)w. 

We have the following identities: 

(3.11) (A 0 ( w + u)9 t w, w) = ^d t (A 0 fw + u)w, w) - ^(d f A 0 (w + u)w, w); 

(3.12) (Ai(w + u)5 x w, w) = ^d x (^i(w + u)w,w) - ^(d x Aifw + u)w, w). 

We consider (13.1011 and take the inner product with w, which yields 

(3.13) (7lo( w + u)<9tw, w) + (Ai(w + u)5 x w, w) = (Di(u, u, w)w, w). 

Using the identities (I3.11H and (13.1211 in (13.13(1 . we obtain 

^<9t(A)(w + u)w, w) + ^d x {Ai{ w + u))w, w) 

(3.14) = ^(d t A 0 (w + u)w, w) + ~ (d x Ai (w + u)w, w) + (Di(u,u,w)w,w). 
Therefore, if we integrate equality (13.131) over R x [0, T], we have 


\ [ (A 0 (w(T) + u)w(T), w(T)) dx - ]- f (A 0 (w(0) + u)w(0), w(0)) dx 
z Jr z Jr 


(3.15) 


dt 


i^A 0 (w + u) + ^d x A\fw + u) + L>i(u,u,w) 


w, w 


dx =: I. 
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To estimate I, we use (13. 1011 in the following form: 


<9(W = — AL(w + u)<9 x w + D( u, u, w)w. 
d t A 0 ( w + u) = A' 0 (w + u)<9 t w, 


(3.16) 

Then 

(3.17) 
and 

(3.18) 9tAo(w + u) = 7 Lq(w + u)(— j4(w + u)(9 x w + D( u, u, w)w) 
Therefore, we have 

rT 


I = - 


1 0 
r T 


dt f ]-((A' 0 A(w+ u)d x w,w),w) dx 
Jr 2 

+ f dt f ^((^(w + u)d x w, w), w) dx 

Jo Jr] 1 

(3.19) + f dt f Di(u, u, w) dx + f dt f -(A' 0 D(u, u)w • w, w) dx. 

Jo Jr Jo Jr 2 

Then 

7 : [ (A)(w(T) + u)w(T), w(T)) dx - [ dt [ (£>i(u, u, w) dx 
2 Jr Jo Jr 

= \j (A)(w(0) + u)w(0),w(0)) dx + [ dt [ ^((^i(w +u)(9 x w, w), w) dx 
2 Jr Jo Jr 1 

(3.20) 

+ [ dt [ i((AgT»(u,u)w,w), w) - f dt [ i((A' 0 A(w + u)O x w),w,w) dx. 
Jo Jr 2 Jo Jr 2 

Since Aq is positive definite, D\ is negative definite and since Aq, A^, A’ 0 D, A' q A 
are bounded in a neighborhood of the equilibrium point u, we have 


|2 dt 


^ C2 n 

< — w 

- 2 


tjII w ( t )IIo + Cl J 0 ll w Wllc 

(°)llo + 7T f dt f |d x w||w| 2 dx + ^ [ dt [ |d x w||w[ 2 dx 

2 Jo Jr 2 j o Jr 

+77 / dt, I |w||w| 2 dx 

2 Jo Jr 


(3-21) < yl|w(0)||g + y||<9 x w||oo^ dt 

where c, c±, c 2 , c 3 , C4, c 5 , c 6 G M + . 


11 2 , c 5 11 I, f ,, || | ,2 

l w Ho + TtIMIoo / dt | |w| | n 
2 Jo 


0 ) 


The embedding of H 1 in L^, where a is the constant of the embedding, yields 
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Thus, from the definition of functional N 2 (t) in (I3.5j) . we have 

(3.23) 11 d x w11oo < a(||5 x w|| 0 + ||<9 xx w|| 0 ) < N 2 (T), 
and 

(3.24) IMIoo < a(||w|| 0 + ||5 x w|| 0 ) < N 2 (T). 

The last term in (13.211) is estimated by 

(3.25) |iV|( 0) + | Nl(T) + | Nl(T). 

So, using (13.341) . we have 

(3.26) |||w(T)||2 + | £ 11 w(t) | |q dt < C (IVf (0) + IVf (T)), 
and, therefore 

(3.27) l|w(T)||o + [ T 11w(t)||g dt < Ci(N$(0) + N 2 (T)). 

Jo 

□ 

Let us now estimate the first and second order derivatives. 

Lemma 3.2. If N 2 (T) < e, Then, for l = 1,2, 

(3.28) sup o <t<T |H*)|If + [ T |Hi)|If dt < C(Af(0) + Af(T)). 

Jo 

Proof. Apply the first space derivative to system (13.41) and take the inner product 
with <9 x w, which provides 

(3.29) (9 x (A 0 (w + u)5 t w +Ai(w + u)«9 x w) • <9 x w = [d x (A 0 D)w + (A 0 D)d x w] ■ d x w. 
We have the following identities: 

3 x (A 0 (w + u)<9 f w) • d x w = ^d t ((A 0 {w + u)d x w) • <9 x w) - ^(9fA 0 5 x w) • <9 x w 

(3.30) + (d x A 0 d t w) ■ d x w: 


(3.31) <9 x (Ai(w + u)3 x w) -d x w = ^9 x ((Ai(w + u)3 x w) • <9 x w) + ^(<9 x Ai<9 x w) • <9 x w. 

If we integrate equality ()3.29l) over R and use the previous identities, the term 

<9 x ((H(w + u)<9 x w) • 3 x w) 

vanishes and then we have 


1 d_ 

2 dt 


(A 0 (w + u)<9 x w) • 8 x mv dx - / ((A 0 L>)<9 x w) ■ d x w dx 


(3.32) = f | ^d t A 0 d x w - d x A 0 d t w - ^<9 x Ai<9 x w 1 • d x w + (d x (A 0 D) w) • <9 x w dx. 
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To estimate the right-end side of (13.321) . we use (13.41) in the following form: 
(3.33) <9 t w = —Ad x w + D'w. 

Then 


dtAo = AIq^w = A' 0 (—Ad x w + Dw). 
Thus, equality (13.321) is 

1 d 

2 dt 


■ / (Ao(w + u)<9 x w) • d x w dx - (( A 0 D)d x w ) • d x w dx 

J M J M 

r i i i 

/ --((A' 0 Ad x w, d x w),d x w) + -(A' 0 d x wAd x w, <9 x w) + -(A’ 0 Dwd x w, d x w) dx 
Jr * 1 1 


1 


(3.34) 

Using (13.341) . we have 


--((A)^9xW,<9 x w),3 x w) + / (A 0 D'(d x w)w,d x w) dx. 

£ ./TU 


N3rw(T)||§+ f \\d x w{t)\\l dt 
Jo 

< C'i(e) ^||9 x w(0)||o + j j |3 x w| 2 |d x w| etedf'j 
(3.35) < C'i(e)M|9 x w(0)|[o + (||w|| 00 +||5 a; w|| 00 )y ||9 x w(i)||o dt j . 


In the same way, we perform the second space derivative of (13.101) and take the 
inner product with d xx w, which provides 


(3.36) d xx (A 0 d t w + A x d x w) ■ d xx w = d xx (A 0 Dw) ■ d xx w. 

We have the following identities: 


(3.37) 


d xx (Aodtw) ■ d xx w 


(3.38) 


= ^d t ((A 0 d xx -w) • d xx w) - ^(d t A 0 d xx w) ■ d xx w 


(3.39) 


+ 2(d x A 0 d xt w ) • d xx w + (d xx A 0 d t w) ■ d xx w ; 


<9xx(-4i(W)<9 x w) • d xx w 










14 


R. BIANCHINP & R. NATALINI* 


If we integrate ()3.36l) over R and we use the previous identities (I3.39I) - (I3.40I) . the 
term 

j x (t -'3 1 (txx) ‘ *9 xx w) 

vanishes. 


In this way, we have 
Id 
2 dt 


(A 0 d xx w) • d xx w dx- (. A 0 Dd xx w , d xx w) dx 

Jr 


7 }d t A 0 d xx w - 2d x A 0 d xt w - d xx A 0 d t w > • d xx w dx 


-d x Aid xx \t+d xx Aid x w > -d xx w+(d xx (AQD)vf, d xx w)+((A 0 D)d x w, d xx w) dx 


- ^(A’qAOx w,<9 xx w) + A' 0 Dw,d xx w ) + (A' 0 d x wA'd x w,d x w) 


+- A' {) d x yj Ad xx ys — 2(^4 / 0 9 x wL) / 9 x w, w) — 2A' G d x wDd x 'w — A! Q d xx vrD-w 


(3.41) - (A 0 A"d x w • d x w. d x w) - (A 0 A'd xx w, d x w) + (A' 0 D'd xx w, w) | • d xx w dx. 
Therefore 

\\d xx w(T)\\l+ [ \\d xx w(t)\\l dt 

Jo 

< C 2 (e)< ||d xx w(0)|||j + [ [ (\d x w\ 2 + |d xx w| 2 )(|d x w| + |w| + |9 x w|) dxdt 

[ Jo J [0,1] 

(3.42) 

< C 2 (e) j||d xx w(0)||o + (llwlloo + ||9 x w||oo) J (||9 x w(t)||o + ||9 xx w(f)||o) dt 

□ 

We note that, by Lemma \3.1\ and Lemma \3.“A we are able to prove inequality 
(13.61) . Then, Proposition^ 7] is proved. 

To conclude the proof of Theorem \3. /[ we are going to prove the statement about 
the asymptotic behavior of solutions. 

Let us set 

(3.43) E (t) = ^||w(t)|||. 

Taking the time derivative of (13.201) and using the embedding of H 2 ( R) in L°°(R), 
we have 
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(3-44) T^IMOIIo + Ci||w(t)||§ < C 2 11w(t)11 2 11w(t)11o• 

Taking the time derivative of (13.341) . we obtain, in the same way, the estimate on 
the first derivative of w. Therefore 

(3-45) ~\\dMt)\\ 2 o + c 3 \\dMm 2 o < c 4 ||w(t)|| 2 ||a r w(t)||2. 

Finally, from the time derivative of (13.411) and using Morrey’s Theorem, we have 
the second order estimate 

(3-46) ~\\d xx w(t)\\l + c 5 \\d xx w(t)\\ 2 0 <c 6 ||w(i)|| 2 ||a xx w(i)||g. 

We note that ci, C 3 , C 5 depend on the norm of Aq, D, while c 2 , C 4 , c$ depend on the 
norm of Aq, A' 0 , A, A', D, D'. Summing (13.441) . (I3.45[) and (|3.46[) . from (I3.43[) we have 

(3.47) d t E + fiE < pE 3/2 , 

where p = c 1 + C 3 + C 5 and v = c 2 + C 4 + cq. 

Now, thanks to the assumptions, we can choose the initial data small enough such 
that 

E (0) < (2) , 

for some 0 < 7 < min{g, ^}. From the smallness of the initial data, we can obtain 
the inequality 

pE 3/2 < 7 E, 

at least for small times. Therefore, using (13.471) we have 

<9 f E < (7 - n) E. 

and <E2D follows directly from Gronwall’s Lemma, taking ft = v — 7 . 
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